SYMMETRY REDUCTIONS OF A NONLINEAR OPTION 
PRICING MODEL 

L. A. BORDAG* 

Abstract. The studied model was suggested to design a perfect hedging strategy for a large 
V^ , trader. In this case the implementation of a hedging strategy affects the price of the underlying 

f^^ 1 security. The feedback-effect leads to a nonlinear version of the Black-Scholes partial differential 

f^ ' equation. Using the Lie group theory we reduce the partial differential equation in special cases to 

^vj ' ordinary differential equations. The found Lie group of the model equation gives rise to invariant 

solutions. Families of exact invariant solutions for special values of parameters are described. 
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("j . 1. Introduction. In a series of works 0) Eli 01 and jH] a model for a liedging 

strategy in an illiquid market was suggested. In the model the implementation of 
(— I . a hedging strategy affects the price of the underlying security. For a large trader 

a hedge-cost of the claim differs from the price of the option. The feedback-effect 
leads to a nonlinear version of the Black-Scholes partial differential equation, 

2 [1- pS\{S)ussr 
> . 

r^ i with S G [0,oo), t e [0, r]. As usual, S denotes here the price of the underlying 

^— ^ ' asset and u{S,t) denotes the hedge-cost of the claim with later defined payoff, 

^ I . which is different from the price of the derivatives product in illiquid markets, t is 

f— ^ ' the time variable, a defines the volatility of the underlying asset, p is a measure for 

\^ . the feedback-effect of a large trader, X{S) is chosen in a way to obtain the desired 

^D ' payoff. The values of p and X{S) might be estimated from the observed option 

prices. In dependence on the propositions on the market different variations of 
the Black-Scholes formula can accrue like in a well known model ^21- Usually the 
volatility term in the Black-Scholes formula will be replaced to fit the behavior of 
the price on the market. The modeling process is not finished now and new models 
can appear. An analytical study of these equations may be useful for an easier 
classification of models created. 

Frey and co-authors studied equation 1)1. i|l under constrictions and did some numer- 
C^ ' ical simulations. Our goal is to investigate this equation using analytical methods. 

We study the model equation (|l.l(l using methods of the Lie group theory. This 
method has a long tradition beginning with the work of S. Lie 0. The applications 
of this method are connected with an obvious limitation of group-theoretical meth- 
ods based on local symmetries because many nonlinear partial differential equations 
do not have local symmetries. The modern description of the method and a large 
number of applications can be found in ^T] , ^U] , ^21 > [Z| ^ [HI - 

In Section [21 we find the Lie algebra and finite equations for the symmetry group 
of equation Hl.l() . For a special form of the function X{S) it is possible to find two 
functionally independent invariants of the symmetry group. 

Using the symmetry group and its invariants we reduce the partial differential equa- 
tion Hl.l|l in special cases to ordinary differential equations in Section |21 We study 
singular points of the reduced equations in Section ^ and describe the behavior of 
invariant solutions. For a fixed set of parameters the complete set of exact invariant 
solutions is given. 
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2 L. A. BORDAG 

2. Lie group symmetries. Let us introduce a two-dimensional space X of 
independent variables {S, i) G X and a one-dimensional space of dependent variables 
u G U. We consider the space C/(i) of the first derivatives of the variable m on S' and t, 
i.e., {us, Ut) G C/(i) and analogously we introduce the space C/(2) of the second order 
derivatives {uss,ust,utt) e C/(2)- We denote hy M = X x U a. base space which 
is a Cartesian product of pairs {x,u) with x = {S,t) (z X, u € U. The studied 
differential equation Hl.l(l is of the second order and to represent this equation as 
an algebraic equation we introduce a second order jet bundle M^^^ of the base space 
M. This space has the form 

M(2) =X xU X t/(i) X C/(2) (2.1) 

and possesses a natural contact structure. We label the coordinates in the space 

M(2) by w = (S, t, u, us, Ut, uss, ust, utt) e M^^) . 

In the space M^^^ equation (|l.ll) is equivalent to the relation 

A(w) = 0, weM^^\ (2.2) 

where we denote by A the following function 



A{S,t,U,Us,Ut,Uss,Ust,Utt) =Ut-\ ::;—— ^c<\rc<\.. \2- (2-3) 



CT^S^ uss 

2 {l-pSX{S)ussy 



We identify the algebraic equation (|2.2(l with its solution manifold La defined by 

La = {w& M^^^\A{w) = 0} C M'-^\ (2.4) 

Let us consider an action of a Lie-point group on our differential equation and its 
solutions. We define a symmetry group Ga of equation (|2.2() by 

GA = {geDiff(M(2))| g: La^La}, (2.5) 

consequently we are interested in a subgroup of Diff (M^^^) which is compatible with 
the structure of La- 

As usual we first find the corresponding symmetry Lie algebra Vif f a{M^'^'') C 
'Diff{M^'^^) and then use the main Lie theorem to obtain Ga and its invariants. 
We denote an element of a Lie-point vector field on M by 

V^aS,t,u)-^+T{S,t,u)^^+cp{S,t,u)^, (2.6) 

where ^{S,t,u),T(S,t,u) and <f>[S,t,u) are smooth functions of their arguments, 

V G Viff{M). 

If the infinitesimal generators of 5 £ Ga exist then they have the structure of the 

type H2.6|l and form an algebra 2?i///^(M). 

A Lie group of transformations acting on the base space M induce as well the 

transformations on M^'^\ 

The corresponding algebra I?i//^(M*^^') will be composed of vectors 

pA^'^v = as, t^^)gg + ^is, ^' ")^ + ^(s, t, U)— 

+ cl)^{S,t,u)j—+cj,\S,t,u)^— (2.7) 

ous out 

ouss oust OUtt 
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which are the second prolongation of vectors V. Here the smooth functions (f)^ {S,t,u) ,| 
(p^{S, t, u), (l)^^{S, t, u), (?!)"^*(5, t, u) and 0"(S', t, u) are uniquely defined by the func- 
tions ^{S,t,u),T{S,t,u) and (f){S,t,u) using the prolongation procedure (see [TT] . 

Cni, El, 0, El). 

Theorem 2.1. The differential equation |J. j|) wii/i an arbitrary function X{S) 
possesses a trivial three dimensional Lie algebra Dif f^(AI) spanned by generators 

d d d 

ot ou au 

Only for the special form of the function X{S) = uiS'^, where uj,k Q R equation ^1.1}) 
admits a nontrivial four dimensional Lie algebra spanned by generators 

d d d d d 

ot ou on Ob on 

Proof. The symmetry algebra Viff^ ( Af '^^ ) of the second order differential equation 
()2.2|) can be found as a solution of the defining equations 

pr^^^V{A) =0 {mod{A^O)), (2.8) 

i.e., the equation 1)2. 8|l should be satisfied on the solution manifold La- 

For our calculations we will use the exact form of the coefficients 0*(5, t, u) and 

(j)^^{S,t,u) only. The coefficient ip*{S,t,u) can be defined by the formula 

(j)*{S, t, u)^(f)t+ Ut(t)u - Usit - UsUtiu - UtTf - {ut)^Tu, (2.9) 

and the coefficient (f)^' ^ {S , t , u) by the expression 

</)^^(5, t, u) = <j)ss + 2us4>su + uss<Pu (2.10) 

+ ("s)^0«u - 2uss^S - UsCsS - 2(u5')^^Su 

- iususs^u - {us)^^uu - 2ustTs - utTss 

- 2usUtTsu - (utuss + 2usust)Tu - {usYutTuu, 

where the subscripts by ^, r, 4> denote corresponding partial derivatives. 
The first equations of the set p. 8(1 imply that if F € Vif f ^{M) then 

S,{S,t,u) ~ aiS, T{S,t,u) — 02, (j){S,t,u) = a^S + Oi + a^u, (2-11) 

where 01,02,03,04,0^ are arbitrary constants and ^,r, (f> are coefficients in the 

expression (|2.6() . 

The remaining equation has a form 

aiSXs{S)-{ai~a5)X{S)^Q. (2.12) 

Because this equation should be satisfied for all S identically we obtain for an 
arbitrary function X{S) 

01=05=0, -^ ^{S,t,u) ^ 0, T{S,t,u)—a2, (l){S,t,u) — a^S + 04. (2.13) 

Finally, T>iff^{M) admits the following generators 

with commutator relations 

[14, 1^2] = [14, ■^^3] = [V2,V3] = 0. (2.15) 
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If the function X{S) has a special form 

A(S')=w5^ uj,kGR (2.16) 

then the equation H2.12|l on the coefficients of (|2.6I) is less restrictive and we obtain 
^{S,t,u) = aiS, T{S,t,u)—a2, (l){S,t,u) ^ {I — k)aiu + a^S + a4. (2.17) 
Now the symmetry algebra T>iff^{M) admits four generators 

with commutator relations 

[1^1, 1^2] = [■^^i, ■^^3] - [^1, 14] = [V2, V3] - 0, 

[V2, V4] = -kV2, [V3,V4] = (1 - k)V3. (2.19) 

D 

Remark 2.1. In the general case the algebra (77) possesses a two dimensional 
Abelian sub-algebra. For the cases fc = 0, 1 the Abelian sub-algebra is three di- 
mensional \2.19\) and we see later that the corresponding equations iV.<S)) became 
autonomous. 

The symmetry algebra Vif f ^{M) defines by the main theorem of S. Lie [5| the 
corresponding symmetry group Ga of the equation (|2.2(l . To find the closed form of 
transformations for the solutions of equation (|l.l|l corresponding to this symmetry 
group we just integrate the system of ordinary differential equations 

^=C(5,t,w), (2.20) 

j^=T{Srt,u), (2.21) 

dii ~ 

— =^{S,t,u), (2.22) 

with initial conditions 

SU=o = S, iU=o = t, uU=Q = u. (2.23) 

Here the variables S, t and u denote values S, t, u after a symmetry transformation. 

The parameter e describes a motion along an orbit of the group. 

Theorem 2.2. The action of the symmetry group Ga of fl.lp with an arbitrary 

function X{S) is given by \2.24^ - \2.2f)]) . If the function \{S) has the special form 

\2.16]) then the symmetry group Ga is represented by i2.2T\ )- l2.2fJ\} . 

Proof. The solutions of the system of ordinary differential equations H2.20|l with 

functions £,,t, (j) defined by (|2.13|) and initial conditions ()2.23f) have the form 

S^S, (2.24) 

i=t + a2e, (2.25) 

u ~ u + a^Se + a^e, ee(— 00,00). (2.26) 

The equations (|2.24l) - (|2.26|l are the finite representation of the symmetry group Ga 
which corresponds to the symmetry algebra defined by H2.14|l in case of an arbitrary 
function X{S). 
If the function X{S) has a special form given by l|2.16|l we obtain a reacher symmetry 
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group. The solution of the system of equations H2.20|) with the functions ^,r, (f) 
defined by H2.17|l and initial conditions (|2.23|l have the form 

S^Se''^", ee (-00,00), (2.27) 

i = t + 026, 

aik 

+ (T .^ (^°'^'''^' - 1)' M 0, A; ^ 1 (2.28) 

ai(l — k) 

u = ue"'" + a-iSee"'' + —{e"^'' - 1), fc = 0, 
ai 

u = u+ —Sie"'" - 1) + 046, fc = 1, (2.29) 

where we assume that ai 7^ because the case ai = coincides with the former 

case (|T^ - (|T^ . n 

We will use the symmetry group Ga to construct invariant solutions of equation 

(|1.1|) . To obtain the invariants of the symmetry group Ga we exclude e from the 

equations (|2.24|) - (|2.2t)|l or in the special case from equations (|2.27|I - H2.29|) . 

In the first case the symmetry group Ga is very poor and we can obtain just the 

following invariants 

invi = S, (2.30) 

inv2 = u — (035 + 04)702, 02 7^ 0. 

These invariants are useless because they do not lead to any reduction of p.l|) . 
In the special case (|2.16f) the symmetry group admits two functionally independent 
invariants of the form 

mwi = log S" + ai, a = ai/a2, a2 7^ (2.31) 

inv2 ^ u S^'^-'^l (2.32) 

In general the form of invariants is not unique because each function of invariants 
is an invariant. But it is possible to obtain just two non trivial functionally inde- 
pendent invariants which we take in the form H2.31|l . (|2.32() . The invariants can 
be used as new independent and dependent variables in order to reduce the partial 
differential equation (|1.1(1 with the special function A(iS') defined by (|2.16|1 to an 
ordinary differential equation. 

3. The special case X{S) = ujS'^. Let us study a special case of equation 
(|1.1() with X{S) ~ LoS^ , k G R. The equation under investigation is now 

, cr^S'^ uss n f'i^\ 

"' + ^- (1 - bs'^+^ussr ^ ^ ^^-'^ 

with the constant b = pto. As usual we suggest that p e (0, 1). The value of the 
constant uj depends on the corresponding option type and in our investigation it 
can be assumed that oj is an arbitrary constant, cj 7^ 0. The variables S, t are in 
the intervals 

5>0, te [0,T], r>0. (3.2) 



Remark 3.1. The case b ~ Q, i.e. p — Q or uj — Q leads to the well known linear 
Black-Scholes model and we will exclude this case from our investigations. 
We will suppose that the denominator in equation Hl.l|l (correspondingly H3.1|l ) is 
non equal to zero identically. 
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Let us study the denominator in the second term of the equation H3.1[l . It will be 
equal to zero if the function u{S,t) satisfies the equation 

1 - bS''+\ss = 0. (3.3) 

The solution of this equation is a function uq(S, t) 

uo{S, t) = —l—S'-'^ + Sc,{t) + C2{t), 5 ^ 0, M 0, 1, 

un{S,t) = -hog S + Sci{t)+C2{t), b^O,k = l, (3.4) 



uo{S,t)^ ]-SlogS + Sci{t)+c2{t), b^Q,k^O, 



where the functions ci(t) and C2{t) are arbitrary functions of the variable t. 

Subsequently we will suggest that the denominator in the second term of the equa- 
tion 13. 1|) is not identically zero, i.e., a solution u{S,t) is not equal to the function 
uo{S, t) l|3.4|l except in a discrete set of points. 
Let us now introduce new invariant variables 

z — log S + at, a 7^ 0, 

v = uS^''-^\ (3.5) 

After this substitution equation (|3.1I) will be reduced to an ordinary differential 
equation 

"^- + T (l-5(... + (l-2fc)..-fc(l-%))^ ^"' '^''^'- ^'-'^ 

Elementary solutions of this equation we obtain if we assume that v = const, or 
Vz = const. . It is easy to prove that there exists the trivial solution w = if A: 7^ 0, 1, 
and the solutions v = const. ^ Q, v — const. 7^ if fc = 0, 1 only. The condition 
that the denominator in (|3.6|l is non equal to zero, i.e., 

(1 - h{v,, + (1 - 2k)v, - k{\ - k)v)f ^ (3.7) 

corresponds to equation (|3.3|l in new variables z,v. 

If the function v{z) satisfies the inequality H3.7|) then we can multiply both terms of 
equation (|3.()|) with the denominator of the second term. In equation (|3.()|) all coeffi- 
cients are constants hence we can reduce the order of the equation. We assume that 
v,Vz ^ const, and choose as a new independent variable v and introduce as a new 
dependent variable x{v) = Vz{z). This variable substitution reduces equation H3.6|l 
to a first order differential equation which is second order polynomial corresponding 
to the function x{v)v. Under assumption Ij^.TI) the set of solutions of equation H3.6|) 
is equivalent to a union of solution sets of the following equations 



a; = , (3.8) 



1 k{l-k)v ^Ja^ {a^ -Sabx) 



= -l + 2k- 



Aab^ x"^ bx x Aab^ x^ 

cr2 1 k{l-k)v ^cr2 (cr2 - g a 6 x) 

Aab"^ x^ bx x Aah^ x"^ 



Equations (|3.8|l are of an autonomous type if the parameter k is equal to fc = 0, 1 
only. We see that these are exactly the cases in which the corresponding Lie-algebra 
(|2.19(l has a three dimensional Abelian sub-algebra. The case fc = was studied 
earlier in ^, [5]. In the next section we will study the case fc = 1. 
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4. The special case A — ujS. If we put A: = 1 in l|2.16[l then equation H3.6|) 
takes the form 

where 9 = f^, a, 6 7^ 0. It is an autonomous equation which possesses a simple 
structure. We will use this structure and introduce a more simple substitution as 
described at the end of the previous section to reduce the order of equation. 
One family of solutions of this equation is very easy to find. We just suppose that 
the value Vz{z) is equal to a constant. The equation H4.1|l admits as a solution the 
value Vz = (— 1± Y^)/6 consequently the corresponding solution u{s, t) of (|3.1(l with 
A = LoS can be represented by the formula 

u{S,t) = — {-l±^){\ogS + at) + c, a>0, (4.2) 

puj 

where c is an arbitrary constant. 

To find other families of solutions we introduce a new dependent variable 

y{z)^vz{z) (4.3) 

and assume that the denominator of the equation (|4.1|l is not equal to zero, i.e. 

z 1 

^'(^) 7^-T+ci e^ + C2, i.e. y{z) ^ -- + c^ e\ (4.4) 



where Ci,C2 are arbitrary constants. 

We multiply both terms of equation 14.111 by the denominator of the second term 

and obtain 

yyl^2[y^ + \y- ^'^yz+ [y' + lv+ {^-^y^V^Q, 5^0. (4.5) 

We denote the left hand side of this equation by F{y^ yz). The equation H4.5|l can 
possess exceptional solutions which are the solutions of a system 

dF{y,yz) n t?( \ n ( a c\ 
7. =0, F{y,yz)^0. 4.6 

OVz 

The first equation in this system defines a discriminant curve which has the form 

y{^) = I- (4.7) 

If this curve is also a solution of the original equation (|4.5() then we obtain an 
exceptional solution. We obtain an exceptional solution if g = 4, i.e. a = cr^/8. It 
has the form 

yi^) = l- (4.8) 

This solution belongs to the family of solutions (|4.1()|l by the specified value of the 

parameter q. In all other cases the equation H4.5|l does not possess any exceptional 

solutions. 

Hence the set of solutions of equation 14.5|l is a union of solution sets of following 

equations 

y - 0, (4.9) 

y^{-l±^)/b, (4.10) 



2,1,. 1 .l^'fl \\^ 



y^ = [y'^iy-^-\l2^[ib-y) -y^y^' ^'-''^ 



2,1 Q.I 

y +by-2^+y¥ 



{Tt-y))l^y^' (4.12) 
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where one of the solutions (|4.1U|) is an exceptional solution (|4.8|) by q — A. We denote 
the right hand side of equations (|4.11|) . (|4.12|) by f{y). The Lipschitz condition for 
equations of the type y^ = f{y) is satisfied in all points where the derivative -rJ- 
exists and is bounded. It is easy to see that this condition will not be satisfied by 

y = 0, 2/ = ^, y^oo. (4.13) 

It means that on the lines H4.13|l the uniqueness of solutions of equations H4.11|l . 
(|4.12l) can be lost. We will study in detail the behavior of solutions in the neighbor- 
hood of lines (|4.13|) . For this purpose we look at the equation (|4.5|l from another 
point of view. If we assume now that z, y, y^ are complex variables and denote 

2/(z) = C, Vz{z)^w, C^weC, (4.14) 

then the equation 1)4. 5|l takes the form 

nC,-) = C-^-2(c^ + lc-^)-+(c^ + ^C + ^)c = o, (4.15) 

where 6 7^ 0. The equation (|4.15|) is an algebraic relation in C^ and defines a plane 
curve in this space. The polynomial F{(^,w) is an irreducible polynomial if at all 
roots Wr{z) of -F'(C, Wr) either the partial derivative i^(;(C, Wr) or -F'iu(C, Wr) are non 
equal to zero. It is easy to prove that the polynomial (|4.15() is irreducible. 
We can treat equation (|4.15|) as an algebraic relation which defines a Riemann 
surface T : F{(, w) ~ oi w — w{() as a compact manifold over the ^-sphere. 
The function w{C) is uniquely analytically extended over the Riemann surface F of 
two sheets over the C— sphere. We find all singular or branch points of w(C) if we 
study the roots of the first coefficient of the polynomial F{(, w), the common roots 
of equations 

F{C,w)^0, F^{C,w)^0, CweCUoo. (4.16) 

and the point C = 00. The set of singular or branch points consists of the points 

Ci=0, C2 = ^, (3 = ^- (4.17) 

As expected we got the same set of points as in real case (|4.13|) by the study of the 
Lipschitz condition but now the behavior of solutions at the points is more visible. 
The points (2 , Ca ^re the branch points at which two sheets of F are glued on. We 
remark that 

^C.)4(.-4) + .-i= + ..., t^ = C-|, (4.18) 



where t is a local parameter in the neighborhood of C2 . For the special value of g = 4 
the value w{C2) is equal to zero. 
At the point C3 = 00 we have 




00, 

where i is a local parameter in the neighborhood of C3. At the point Ci = the 
function w{() has the following behavior 

w(C) r^ — ^ — , C -^ Ci = 0, on the principal sheet, (4.19) 

w{C) -- (1 - g) C, C -^ Ci = 0, g 7^ 1, on the second sheet, (4.20) 

w{C) 2b'^C^, C -> Ci = 0, g = 1, on the second sheet. (4.21) 
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Any solution w{Q of an irreducible algebraic equation (|4.15|) is meroniorphic on 
this compact Riemann surface F of the genus and has a pole of the order one 
correspondingly H4.19|l over the point Ci = and the pole of the second order over 
(^3 = oo. It means also that the meroniorphic function w{Q cannot be defined on a 
manifold of less than 2 sheets over the C sphere. 

To solve differential equations H4.11(l and H4.12|l from this point of view it is equiv- 
alent to integrate on F a differential of the type ■:^^^ and then to solve an Abel's 
inverse problem of degenerated type 

^ z + const. (4.22) 



^(0 

The integration can be done very easily because we can introduce a uniformizing 
parameter on the Riemann surface F and represent the integral H4.22|l in terms of 
rational functions merged possibly with logarithmic terms. 

To realize this program we introduce a new variable (our uniformizing parameter 
p) in the way 

C-'-^, (4.23) 

^^(l-rt(.(l+p)-^4)^ 

46(p +1) ^ ' 

Then the equations H4.11(l and H4.12|l will take the form 

The integration procedure of equation (|4.25|l gives rise to the following relations 

2q log (p - 1) + (g - Vg - 2) log ((p + 1) V9 - 2) (4.27) 

+ iq + ^^ 2) log Up + 1)^ +2) = 2{q-l)z + c, q^l,q>0 

2v/F9yarctan ((p + 1)^(^/2) - 2qlog (p - 1) 

+ (2-g)log(4-q(p+l)2) = 2(l-g)z + c, q<0, (4.29) 

where c is an arbitrary constant. The equation H4.26|l leads to 

2glog(p+l) + (g + V9-2)log((p-l)V9-2) (4.30) 

+ (g-^-2)log((p-l)Vg + 2) = 2(g-l)z + c, q^l,q>0 

^ + ilog^^^=z + c, q=l. (4.31) 

p + 1 4 (p + 1)^ 

-2yi^arctan((p - 1)71^/2) - 2glog(l +p) 

+ {2 - q) log {A -q{p- if ) = 2{l-q)z + c, q<0. (4.32) 

where c is an arbitrary constant. 

The relations (|4.27() - (|4.32() are first order ordinary differential equations because of 

the substitutions (|4.14|) and H4.3(l we have 

Ah 

(4.33) 
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All these results can be collected to the following theorem. 

Theorem 4.1. The equation \4-l\l for arbitrary values of the parameters q,b ^ 

can be reduced to the set of first order differential equations which consists of the 

equations 

V, = 0, v, = (-1 ± Vg)/& (4.34) 

and equations \4-.2'J^ - \4-.'J'/J^ . The complete set of solutions of the equation \4.1\j 

coincides with the union of solutions of these equations. 

To solve equations H4.27|l - H4.32|l exactly we should first invert these formulas in 

order to obtain an exact representation p as a function of 2;. If an exact formula 

for the function p — p{z) is found we can use the substitution (|4.33(l to obtain an 

explicit ordinary differential equation of the type Vz{z) — f{z) or another suitable 

type and if it possible then to integrate the final equation. 

But even on the first step we would not be able to do this for an arbitrary value of 

the parameter q. It means we have just implicit representations for the solutions 

of the equation (|4.1|l as solutions of the implicit first order differential equations 

(I07ll - (I0^ . 

4.1. Exact invariant solutions in case of a fixed relation between vari- 
ables S and t. For a special value of the parameter q we can invert the equations 
(|4.27(l and (|4.3()(l . Let us take q — A, i.e., the relation between variables S, t is fixed 
in the form 

2 

z = \ogS+^t. (4.35) 

8 

In this case the equation H4.27|l takes the form 

(p - 1)2 (p + 2) == c exp (3z/2) (4.36) 

and correspondingly the equation (|4.3U|) the form 

(p + lf{p - 2) = c exp (3z/2), (4.37) 



where c is an arbitrary constant. It is easy to see that the equations (|4.36|l and 
(|4.37l) are connected by a transformation 

p^ -p, c ^ -c. (4.38) 

This symmetry arises from the symmetry of the underlining Riemann surface F 
(|4.15() and corresponds to a change of the sheets on F. 
Theorem 4.2. The second order differential equation 

0, (4.39) 



{1 ~ b{v,, - v,)y 



is exactly integrable for an arbitrary value of the parameter b. The complete set 
of solutions for b ^ is given by the union of solutions ^4.4^ , i4-45{ l - i4-4<^ f^^d 
solutions 

3 1 

v{z) ^ d, v{z) = --Z + d, v{z)^-z + d, (4.40) 

where d is an arbitrary constant. The last solution in {4-41^ corresponds to the 
exceptional solution of equation j^.5| ). 

For b — equation i4.39^ is linear and its solutions are given by v{z) — di + 
(i2 exp (3z/4), where di,d2 are arbitrary constants. 
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Proof. Because of the symmetry ()4.38f) it is sufHcient to study eithr the equations 

(I4.36|) or H4.37|l for c G i? or both these equations for c > 0. The value c = can be 

excluded because it complies with the constant value of p{z) and correspondingly 

constant value of Vz{z), but all such cases are studied before and the solutions are 

given by (|4.40() . 

We will study equation H4.37|) in case c ^ R\ {0} and obtain on this way the 

complete class of exact solutions for equations (|4.3t)|) - (l4.37|) . 

Equation (|4.37l) for c > has a one real root only. It leads to an ordinary differential 

equation of the form 

1 2i 

v4^) = -l 



b(2 + ce^ + V4c 



2 + ce 2 -\- ^4 c e "2" -\- c^ e^ ' 

, OO. (4.41) 



b2i 

Equation 14.41|l can be exactly integrated if we use an Euler substitution and in- 
troduce a new independent variable 



r = 2 + ce~ + Y4ce~ -l-c^e^^. (4.42) 

The corresponding solution is given by 



2# 

V{Z) : 



bh + ce^ + ■\/4ce^ +0^3^ A 




2t 



v4 ce 2 + c^ e^^ 



(2 + ce^ + VTce^TcV^j 
62i 



ce 2' -\- ^/4 ce 2 + c^ e^^j 

2i 



(4.43) 



where d e i? is an arbitrary constant. 

If in the right hand side of equation (|4.37|) the parameter c satisfies the inequality 

c < and the variable z chosen in the region 

ze (-00,1 In A) (4.44) 

then the equation on p possesses maximal three real roots. 

These three roots of cubic equation (|4.37() give rise to three differential equations 

of the type w^ = (1 — p^[z))/b. The equations can be exactly solved and we find 

correspondingly three solutions Vi{z), z = 1, 2, 3. 

The first solution is given by the expression 

, . z 2 (2 ( \c\ 3.\\ 

i;i(z) = - - - cos I - arccos II — TT '^ '^ ) ) (4.45) 

- ^ log f 1 + 2 cos f - arccos \\-^-e^ 

16 , / /I / \c\ 3^ 

- g^ log I sm I - arccos I 1 " y ^ 2 
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where d G i? is an arbitrary constant. The second solution is given by the formula 

/ X 2 2 /2 2 / \c\ 3.\\ 

^2(2;) = T - T cos (0^+0 arccos ( -1 + y e 2 II (4.46) 

4 , A „ /I 1 / , |c| 3, 

" 3^ log ( 1 + 2 cos I -TT + - arccos ( -1 + — e ^ 

16 , / /I 1 f , \c\ g 

-- log (^sin(^-. + - arccos (^-l + -e 

where d € i? is an arbitrary constant. The first and second solutions are defined up 
to the point z = | In -jIt where they coincide (see Fig. I4.1f) . 

The third solution for z < jt In -rr is given by the formula 

V3a{z) = ^ - -^ cos f - arccos ( "1 + y ^^ j j (4.47) 

4 , / /I / |c| .' 

— — log 1—1 + 2 cos I - arccos I — 1 + — e 

16 , / /I / , |c| 3. 

~ 36 v ^^^ \ 6 ^'^^'^^^ I ^ y ^ ' 

where d e i? is an arbitrary constant. In case z > | In ^ the polynomial (|4.37|) has 
a one real root and the corresponding solution can be represented by the formula 

vzaiz) = r ~ T cosh i -arccosh ( ~1 + y ^^ j j (4.48) 



log I cosh I -arccosh — 1 + 



3& °V V6 V 2 



TTT log 


-1 + 2 cosh 


—arccosh 


36 ^ V 




yi 



■ e 



c 3. 

-1 + U,^ 



The third solution is represented by formulas v^^2{z) and ^'3.1(2;) for different values 

of the variable z. D 

One of the sets of solutions (|4.43|l . (|4.45|l - (|4.48|l for fixed parameters 6, c, d is 

represented in Fig. 14.11 The first solution (|4.4HI and the third solution given by 

both (|4.47|l and (|4.48l) are defined for any values of z. The solutions vi (z) and V2{z) 

cannot be continued after the point 2 = ^ In -m- where they coincide. 

2 
If we put in mind that z = logS* + ^t and v{z) ~ u{S,t) we can represent exact 

invariant solution of equation (|3.1|) . The solution 14.43|) gives rise to an invariant 

solution u(S,t) in the form 



u{S, t) = -— ( 1 + cSh^* + \/2 c5ie^* + c^ S^e- 
ujp 



-— fl + c^^e^* + V2c^«e^* + c2^3e^* 1 (4.49) 
ujp\ 



-Aiog( ( l + cS'5eT^« + y2cS'te^* + c2S'3eT* 



1 + cS'Se-nr* + V2cS'§e^* + 02 536^-' 
where d e i?, c > 0. 
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Fig. 4.1. Plot of the solution v{z), \4.4'^ , (thick solid line), v\{z), 14-45^ , (short 
dashed line), V2{z), |4.^6| ), (long dashed line) and the third solution V3^i(z),V3^2(z), 14-4T\I > 
i4-4'^ > which is represented by the thin solid line. The parameters takes the values \c\ = 
l,q = A,d = 0,b = 1 and the variable z 6 (—5, 4.5). 



In case c < we can obtain correspondingly three solutions if 

0<5< ( ^ ) expf -^t 



(4.50) 



The first solution is represented by 



Mi(5, t) = ^ ( logS'+ ^— t ) ^ cos ( - arccos 1 ~ — 5*2 6 is * 

ujp \ 8 J Lup \3 \ 2 

4 , /, „ /I /, \C\ „3 3o£, 

log 1 + 2 cos - arccos I 1 — ^-^ S^e ^« 

Sujp \ \3 \ 2 

16 , / . A f, |c| „3 3^; 
log sm - arccos 1 — ^-^ S^e i<i 



(4.51) 



where d d R, c < 0. The second solution is given by the formula 

^2 



U2{S,t) = — flog 5+ ^t 
ojp 



2 /2 2 / , |c| ^3 3^ 

— cos -TT H — arccos — IH b^e le 

ujp \6 3 \ 2 



4 , A „ /I 1 / \c\ „3 3.% 

- — log 1 + 2 cos -vr + - arccos -1 + V S'^e^^* 
3ojp \ \3 3 \ 2 



(4.52) 



16 



1 1 



log sm -TT + - arccos [ ~1 + '-^ S^e ^f^ 

3ujp \ \6 6 \ 2 

where d e i?, c < 0. The first and second solutions are defined for the variables 
under conditions H4.50|) . They coincide along the curve 



and cannot be continued further. 
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Fig. 4.2. Plot of solutions u{S,t),ui(S,t),U2{S,t),U3,i{S,t),U3^2{S,t) for the param- 
eters \c\ = 0.5, q = 4, 6 = 1.0, d = 0. The variables S,t lie in intervals S G (0,9) and 
t g [0,2.0]. All invariant solutions change slowly in t-direction. 



The third solution is defined by 

U3.i{S,t) = — logS'+— t cos - arccos -1 + V -S^'e-nr* 

ujp \ 8 J ivp \S \ 2 

_ ^iog('-l + 2cos('i arccosf'-l + M^le^A')') (4.53) 

16 , / /I / _, \C\ „3 3^, 

log cos - arccos l-l + ^-^S^e^o 

Sup \ \6 \ 2 

where d £ R and S", i satisfied the condition (|4.5U|I . 

In case log S + ^t > | In ^ the third solution can be represented by the formula 

U3 2(S,t) = — (log 5*+ ^t] coshf -arccoshf -1 + i^S'leT^r* 

ojp \ 8 J Lup \3 \ 2 

^^ log f cosh fiarccoshf-l + M^le^A)) (4.54) 



3cup \ \6 V 2 

log (-1 + 2 cosh ( -arccoshf -l + ^S'5e^* ] ] ] + d. 



3ujp \ \3 \ 2 

The solution u{S,t) H4.49|l and the third solution given by U3_i, -^3,2 H4. 5311 . 14. 541) 
are defined for all values of variables t and S > 0. They have a common intersection 
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curve of the type S = const. exp(— ^t). The typical behavior of aU these invariant 

solutions is represented in Fig. 14.21 

Previous results can be summed up in the following theorem describing the set of 

invariant solutions of equation 1)1. 1|1 . 

Theorem 4.3. 

1. The equation ^.1}) possesses invariant solutions for the special form of the 
function X{S) given by \2.1(^) only. 

2. In case ^2.16]) the invariant solutions of equation ^1.1\) are defined by or- 
dinary differential equations 1^3. 8\l . In special cases fc = 0, 1 equations ^3.8\) are of 
an autonomous type. 

3. If X{S) = ujS, i.e. fc — 1, then the invariant solutions of equation U.l\l 
can be defined by the set of first order ordinary differential equations \4.2^/\ l~ \4.32^ 
and equation ^.34}^ . 

If additionally the parameter q — A, or equivalent in the first invariant ^2.31\l we 
chose a — cr^/S then the complete set of invariant solutions U.l\l can be found 
exactly. This set of invariant solutions is given by formulas {4.4^ - ^4-54)) o.'nd by 
solutions 

u{S,t)^d, u{S,t)^~3/b{\ogS + ah/8), u{S,t) ^ l/b {log S + cr^t/8), 

where d is an arbitrary constant. This set of invariant solutions is unique up to the 
transformations of the symmetry group Ga given by theorem \2.Sl 
The s olutions u{S, t) gi^, ui(5,i) ^^,U2{S,t) (|i33), wg.iCS*,*) ^^J^ .u-,o(S.t) 
()4.54(l . have no one counterpart in a linear case. If the parameter p —> then 
equation Hl.l|) and correspondingly equation (|3.1|l will be reduced to the linear 
Black-Scholes equation but solutions H4.49|) - (|4.54|) which we obtained here will be 
completely blown up by p ^ because of the factor 1/& = \/{ujp) in the formulas 
(I4.49(l - (|4.54(l . This phenomena was described as well in j^, for the invariant 
solutions of equation 1)3.1(1 with fc = 0. 
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